Zoological Journal of the Linnean Society

The Dinaric karst of Croatia encompasses a network of over 10,000 caves and represents one of the worlds most important subterranean biodiversity hotspots. It is inhabited by remarkably diverse and often endemic species, including planarian flatworms, which are among the rarest macroinvertebrates encountered in cave habitats. Although the presence of cave planarians has long been known, no integrative research on this group has been conducted to date, and the evolutionary relationships between these animals and their surface water counterparts are currently unresolved. To address these gaps, we combined field sampling, phylogenetic analysis based on COI and 18S genes, and phenotypic characterization. Our results show that cave planariids in Croatia belong to at least three genera and are more widespread and diverse across both Croatia, and the broader Dinaric karst, than previously assumed. We increased the number of cave records in the Dinaric karst from 26 to 37 and documented cf. Atrioplanaria and Phagocata in Croatian caves for the first time. Phylogenetic reconstructions suggest numerous independent cave colonization events, including multiple instances within the genera Crenobia and cf. Atrioplanaria. Variation in pigmentation and eye reduction, both within and between populations, further reveal heterogeneous evolutionary trajectories of cave-associated phenotypes. The biogeographical patterns and high genetic diversity we report here point to a complex evolutionary history of planariids in the Dinarides. Our newly generated molecular phylogenies and systematic documentation of trait variability establish Planariidae as a valuable model for studying mechanisms underlying convergent evolution of pigment loss and eye reduction in cave environments.

Color polymorphism can facilitate local adaptation, maintain intraspecific diversity, or reflect early stages of speciation. Clownfishes (Amphiprion spp.) typically display a simple black, orange, and white pattern, but the saddleback clownfish (Amphiprion polymnus) shows striking variation in melanism and the number of vertical bars, which are thought to play a role in species recognition. In 2024, a revision on iNaturalist split A. polymnus into multiple species based solely on color pattern and geographic range. This raises the question of whether these morphs represent true species or intraspecific polymorphism, which we tested using genomic and image-based data. We sampled 97 individuals from seven populations across the species range and quantified color patterns from standardized photographs. Phenotypic and genomic analyses reveal a complex pattern of divergence. Image analysis identified three distinct phenotypic clusters, with A. polymnus, A. annamensis, and A. laticlavius each showing consistent differences in saddle shape and vertical bar extent. ADMIXTURE resolved three distinct genetic groups corresponding to the morphs. Pairwise FST (0.54-0.71) and dxy indicate extremely high differentiation between A. polymnus and A. annamensis, consistent with species-level divergence, whereas A. laticlavius shows much lower differentiation from A. polymnus (FST 0.09-0.18) and higher differentiation from A. annamensis (FST 0.64-0.66). Overall, phenotypic and genomic data show structured variation, but the status of A. laticlavius remains ambiguous. Our study reveals clear and structured divergence across the full range, yet the taxonomic interpretation of this variation remains inherently challenging. The key question remains: do these patterns reflect a single polymorphic species or a complex of closely related species?

The genus Trichiurus is the most economically valuable fish in the family Trichiuridae, currently recognized to include 10 valid species. However, historically numerous morphologically similar congeners have been erroneously assigned as synonyms or subspecies of T. lepturus. In this study, we examined 16 hairtail specimens collected from the southern waters of Java Island, Indonesia. Integrated morphological and mitochondrial phylogenetic analyses (COX1 and 16S rRNA), compared against global Trichiurus sequences, revealed that these specimens form an independent lineage that diverged early from other congeners. Consequently, we describe this lineage as a previously undescribed cryptic species. Diagnostic characters include: first anal-fin spine below 36th-37th dorsal-fin rays; anus below 35th-36th dorsal-fin rays; anteriormost tip of supraoccipital well posterior to posterior distal margin of eye; anterior margin of the pectoral-fin spine non-serrated; fangs on both jaws with barb-like processes; upper jaw long, mean 16.6% (15.5-17.6%) of preanal length; snout short, 12.0% (10.9-13.1%) of preanal length; eye small, diameter 5.3% (4.3-5.7%) of preanal length; and absence of hyperostosis on dorsal cranium. We herein propose the name Trichiurus javaensis sp. nov., and provide a formal morphological description and diagnostic characterization of this species.

The Clouded Tiger Cat Leopardus pardinoides is a recently recognized Neotropical species for which ecological and natural history data are sparse. Knowledge of species distribution and elevational range are largely based upon camera trap studies, and its diet has not been examined. The objective of this study was to better define the geographical and elevational distribution of the subspecies L. pardinoides oncilla in Central America using genetically confirmed records. We also provide the first diet analysis for L. pardinoides. We conducted extensive surveys for scat samples by visual means and with a scent detection dog across mountain ranges in Panama and Costa Rica, confirmed species identity using Sanger sequencing of mitochondrial DNA, and analyzed diet composition via DNA metabarcoding. Altogether, we collected 195 confirmed L. pardinoides scats. Records were strongly associated with high elevations, with a median elevation of 2805 m, and a maximum of 3422 m. Our findings highlight core habitat for L. pardinoides in the Cordillera Talamanca mountains in the border region of Panama and Costa Rica, and we identified isolated populations from the Central Volcanic Cordillera of Costa Rica and the Central Cordillera of Panama. DNA metabarcoding detected 59 vertebrate taxa. Diet was dominated by small mammals (74% of samples), especially cricetid rodents and shrews. Birds were also commonly detected (35% of samples), while reptiles (13%) and amphibians (1.9%) were less frequent. Estimated median adult prey mass was low (25 g), indicating that the species specializes on small prey. Most vertebrate prey identified were species endemic to montane regions of Panama and Costa Rica. These findings provide valuable ecological and natural history data on Clouded Tiger Cats in Central America that will inform conservation and management of this recently recognized species.

The diversity of body shapes is one of the most prominent features of phenotypic variation in mammals. Yet, mammalian body shapes are poorly quantified and the underlying components contributing to its diversity as well as its relationship to other components of the skeleton are rarely tested. Here, we use lagomorphs (hares, rabbits and pikas) as a model system to (1) investigate which components of the skeleton contributed the most to body shape diversity, (2) examine the relationships between body shape and relative limb lengths, and (3) test how body size, ecotype, burrowing behavior, and locomotor mode influenced variation in lagomorph body shape and appendicular morphology. We quantified the body shape and functional proxies of the appendicular skeleton in 40 lagomorph species from osteological specimens held at museum collections. Using phylogenetic comparative methods, we found the relative length of the ribs and elongation or shortening of the thoracic and lumbar regions contributed the most to body shape evolution across lagomorphs. Second, we found that only leporids (hares and rabbits) exhibited a significant relationship between limb length and body shape, where more elongate species exhibit relatively shorter forelimbs and hindlimbs. Lastly, we found that models incorporating body size were the best predictors of lagomorph body shape and the majority of the appendicular traits, whereas models incorporating burrowing behavior and locomotor mode were largely poor fits. Broadly, these results indicate that larger lagomorphs tend to exhibit more robust body shapes with longer, more gracile forelimbs, whereas smaller lagomorphs tend to exhibit more elongate body shapes with shorter, more robust forelimbs. Overall, this work contributes to the growing understanding of mammalian body shape evolution and demonstrates the importance of not omitting body size in ecomorphological analyses.

Our study provides the first comprehensive karyotype evolution analysis of the insect order Zoraptera. We present karyotypic descriptions of seven species across two families: Zorotypidae (Usazoros hubbardi and two Zorotypus spp.) and Spiralizoridae (Centrozoros gurneyi, Spiralizoros magnicaudelli, and two Spiralizorose spp.). These results facilitate a critical evaluation of existing cytogenetic knowledge in Zoraptera and the evolution of karyotypic traits across Polyneoptera. Most notably, we refute the presence of holocentric chromosomes in Zoraptera. Also, we show that the XY sex chromosome system is prevalent and likely ancestral within the order. Furthermore, by integrating the chromosome numbers of the studied species with a dated molecular phylogeny of Zoraptera, we provide the first estimation of the mode of chromosome number evolution for this group. Finally, standard karyotypic features (2n, chromosome morphology, and size) and the distribution of 18S rDNA and (TTAGG) telomeric motif clusters--detected by fluorescence in situ hybridization--reveal highly differentiated karyotypes and genomic structures. This genetic diversity contrasts sharply with the recognized morphological uniformity of Zoraptera.

BackgroundMaximum running speed is a central performance trait, linking morphology, physiology and behaviour to fitness. It is shaped by physical capacity and ecological selection but may also be constrained by ancestry. To examine how these forces interact across macroevolutionary timescales, we conducted an allometric study in a hyper-diverse arthropod taxon--spiders (Araneae). ResultsDrawing on running performance data for 258 species from 64 of the 139 extant spider families, we integrated phylogenetic comparative methods and biomechanical modelling to disentangle the effects of body size, ancestry, leg morphology, ecological guild and preferred locomotor orientation. Maximum running speed varied substantially, both across body mass and among species of similar body mass. By accounting for body mass with a recent biomechanical model, we show that size-specific performance carries a strong phylogenetic signal, and that high-performing runners first evolved within the derived infraorder Araneomorphae.Strong running performance, after accounting for both body size and shared ancestry, was associated with relatively longer legs and, to a lesser extent, ecological guild, but not with leg slenderness or a preference for inverted versus upright locomotion. ConclusionsMacroevolutionary patterns of running performance thus reflect not only variation in body size, but also size-specific leg morphology, ecological differentiation and phylogenetic history. We hope this study contributes to the development of formal evolutionary biomechanics--one that seeks to explain patterns of diversity through the explicit integration of large-scale comparative data, natural history and quantitative models derived from first principles.

This study presents a systematic revision of the suborder Astrophorina for the temperate Pacific coast of the United States and Canada. Major findings include a reduction in the number of species previously thought to range into the region from Japan; validation of most Geodia species erected by Lendenfeld (1910), which were later synonymized by de Laubenfels (1932); the formal description of 10 new species (Poecillastra alaskensis sp. nov., Vulcanella explorata sp. nov., Vulcanella rupta sp. nov., Stelletta cardenasi sp. nov., Stelletta nicolenya sp. nov., Stelletta limuwensis sp. nov., Dercitus (Stoeba) giveni sp. nov., Penares anyapax sp. nov., Penares foxi sp. nov., and Thenea diastra sp. nov.); and one new combination, Penares orientalis comb. nov. Extensive SCUBA-based collection efforts yielded new samples for 11 of the 26 species identified in the region, which enabled an integrative taxonomic approach that combined field photography, fresh material for DNA sequencing, and improved characterization of species ranges and morphological variability in previously described taxa. Illumina sequencing generated complete nuclear ribosomal haplotypes for five species, while Sanger sequencing of the 28S and cox1 loci placed 20 of the 26 species within molecular phylogenies. The use of very short "mini-barcode" amplicons also enabled sequence recovery from historic type specimens up to 137 years old. This study additionally reports the discovery of sponge grounds of abundant, large Geodia at diving depths in Southern California. Together, these results substantially advance our understanding of global astrophorid diversity and systematics, and the biogeography of sponge diversity in the Northeast Pacific. Note about species names: this pre-print is not intended to be a publication of the associated species names for the purposes of zoological nomenclature.

Flight feather vanes are the primary aerodynamic surface of the avian wing. Because loading varies across the wing, vane macrostructure should co-vary with local mechanical demands, yet comparative data on how barb and barbule traits change among remiges and between vane surfaces remain scarce. We quantified barb density, barbule density, barb angle, barb length, and vane width on both vanes at three measurement positions along the rachis of all remiges in four species with contrasting flight modes (white stork, common buzzard, house sparrow, pygmy cormorant), generating over 40,000 measurements across 15 response variables from 992 feathers of 41 individuals. Two complementary generalised additive models characterised variation along the spanwise, inter-vane, and longitudinal axes, and compared outer primaries, inner primaries, and secondaries as functional wing regions. Feather macrostructure varied along all three axes and outer primaries represent the most distinctive region, with lower leading-vane barb density, reduced barb angles, and vane width asymmetry two to three times higher than in inner primaries or secondaries. House sparrow exhibited the densest vane architecture and the highest vane width asymmetry, whereas the low wing-beat frequency species showed complex nonlinear spanwise patterns undetectable by single-feather sampling. Pygmy cormorant barbule density was 39-53% lower than in all other species, matching its wettable plumage strategy. Longitudinal gradients in barb density and barb angle (22-31% decline) were conserved across species. The avian wing is thus functionally regionalised at the macrostructural level, with vane architecture reflecting both aerodynamic and ecological pressures. Summary statementFine-scale vane measurements across all remiges in four species show macrostructural regionalisation of the avian wing, with outer primaries showing the most distinctive vane architecture.

Gorgonopsian therapsids represent a transitional condition in the evolution of synapsid locomotion and postcranial structure. Most descriptions of gorgonopsians have focused on cranial material, however, limiting their usefulness for informing patterns of postcranial evolution on the mammal stem. While some recent work has begun to focus on postcrania, especially the pectoral girdle and forelimbs, comparatively little data are available on the pelvic girdle, hind limbs and tail. We report a new specimen of the late Permian gorgonopsian Aelurognathus tigriceps comprising a partial skull and well-preserved postcranial skeleton, including the near-complete series of dorsal vertebrae and ribs, complete pelvic girdle, hind limbs, feet, and a nearly complete tail. The tail is longer than any other published gorgonopsian. The new material presented here provides an opportunity to better establish broader patterns of morphology in the gorgonopsian postcranial skeleton.

This investigation aimed at compiling all phylogenetic lineages within and around genus Cyanoboletus. The evolutionary inference obtained from the nuclear ribosomal genes internal transcribed spacer region (ITS) suggests that part of the species currently classified in Cyanoboletus belong in lineages separate from the genus, thus suggesting a narrower boundary that includes only the species that develop a strong staining reaction to touch and to air exposure of the context. The separate lineages are the monotypic Cupreoboletus genus and a few species that do not develop such reaction, which are part of a clade together with genera Cacaoporus and Acyanoboletus, thus broadening the concept of Cacaoporus to encompass all of them. The emerging 3C perspective of Cupreoboletus, Cacaoporus and Cyanoboletus offers a remarkably consistent morphological diagnosis, overcoming the problems of a too broad concept for Cyanoboletus. This work reveals that Boletus neotropicus, B. novae-zelandiae and B. sensibilis belong respectively in Cyanoboletus, Cacaoporus and Lanmaoa, and by studying multigene alignment concatenates it identifies lineages that probably represent undescribed species: at least four in Cacaoporus and at least five in Cyanoboletus. Diagnostic tables and dichotomic keys are presented by geographic region. The present work also includes a study of the phylogenetic position of Neoboletus flavosanguineus, a species once classified in Cyanoboletus. The complexity of assigning species epithets in some lineages is addressed, namely for the boundaries between Cacaoporus instabilis and Ca. fagaceophilus as well as the diversity under the names Cyanoboletus sinopulverulentus and Cy. pulverulentus. The overall picture of evolutionary lineages sets a framework for the choice of reference data that can provide, in future phylogenetic studies that involve the 3C, a balanced and efficient coverage. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=197 SRC="FIGDIR/small/724631v1_ufig1.gif" ALT="Figure 1"> View larger version (23K): org.highwire.dtl.DTLVardef@7f618corg.highwire.dtl.DTLVardef@dd6a14org.highwire.dtl.DTLVardef@5f7399org.highwire.dtl.DTLVardef@9e7443_HPS_FORMAT_FIGEXP M_FIG C_FIG

The moonfish Mene maculata (Bloch and Schneider, 1801), the sole extant species of the family Menidae, is widely distributed in the Indo-West Pacific. Although its spawning season has previously been inferred from seasonal variation in the gonadosomatic index (GSI), the histological basis for interpreting GSI as an indicator of spawning-capable status remains limited. Here we describe the gonadal structure of male and female M. maculata and evaluate how GSI and standard length relate to spawning-capable status based on germ cell development. The testis was lobular and exhibited an unrestricted spermatogonial distribution, a structure widely observed among neoteleosts, indicating that the phylogenetic distinctiveness of M. maculata was not associated with a distinctive testicular structure. In the spawning-capable female, oocytes at multiple developmental stages, from primary growth to oocyte maturation, co-occurred within the ovary, indicating asynchronous ovarian development. This finding suggests that M. maculata may be a batch spawner rather than a total spawner as previously inferred. The spawning-capable female had a GSI consistent with previously inferred spawning season estimates. By contrast, histologically examined males were spawning capable at GSI values lower than those previously associated with the inferred spawning season, suggesting that male spawning-capable status may persist beyond the period inferred from elevated GSI alone. This study provides the first histological description of reproductive biology in the phylogenetically distinctive M. maculata and establishes a histological basis for interpreting GSI as a reproductive indicator.

Animal weapons are ecologically important traits that mediate contests over limiting resources and can strongly influence survival and reproduction. Weapon traits often exhibit substantial intraspecific morphological diversity, raising questions about the ecological drivers of this variation. Acrorhagi are weapons produced by sea anemones that are used in intraspecific territorial encounters. Although acrorhagial morphology varies widely within species, patterns of intraspecific variation remain poorly characterized, and the extent to which such variation reflects differences in local intraspecific competition is unclear. Here, we conduct morphometric analyses to characterize within-population variation and allometry in acrorhagial traits of the solitary anemone Anthopleura sola. We show that these traits covary with habitats differing in conspecific density. The number of acrorhagi scaled positively with body size, and individuals occupying a high-density habitat tended to possess more acrorhagi than did similar sized individuals from a low-density habitat. In addition, anemones from high-density habitats exhibited longer acrorhagial cnidae, a pattern that was not explained by differences in body size or acrorhagial density. Together, these results suggest that competitive context influences weapon-related traits at multiple levels of biological organization, potentially via phenotypic plasticity or selective processes. More broadly, our findings highlight how fine-scale ecological variation may contribute to the maintenance of trait diversity within and across species.

Silk glands have been found in two groups of amphipods: the Corophiida and the Ampeliscidae. The silk glands in Ampeliscidae, however, have yet to be examined in detail. Here we report, for the first time, the morphology and distribution of pereopodal glands in the Ampeliscidae, in non-thread producing Synopiidae, and in the Paragammaropsidae. In the Ampeliscidae we found two gland types distributed throughout all pereopods which have the ability to create threads. Pereopods three and four have additional silk extrusion morphology at the tip of the dactylus in which silk is transformed into semi-cylindrical threads used for building domiciles. Synopiid outgroup species have one of the gland types but lack silk extrusion morphology. Using ancestral state reconstruction analysis, we find that glands in the Synopiidae are likely ancestral and hypothesize that silk glands in Ampeliscidae are derived from these ancestral glands. Silk-spinning pereopods in the Paragammaropsidae had similarities with both Corophiida and Ampeliscidae but had distinctions. Ampeliscidae silk-spinning systems bear surprising resemblance to the Corophiida which presents one to reconsider the taxonomic placement of Ampeliscidae and the origins of silk-spinning in amphipods. This is the first comprehensive study on the glandular systems of Ampeliscidae, Synopiidae, and Paragammaropsidae using advanced microscopy, providing pertinent morphological data to the study of arthropod silk gland evolution and complex traits.

Chaetognathifera, a superphylum comprising Syndermata (Rotifera including Acanthocephala), Micrognathozoa, Gnathostomulida and Chaetognatha, is a complex grouping generally recovered as the sister to all other Lophotrochozoa. However, phylogenetic relationships within this group are controversial, in part due to poor sampling, resulting in two key questions. The first is whether Gnathostomulida or Chaetognatha represent the sister group to Syndermata+Micrognathozoa. The second is the phylogenetic position of the former phylum Acanthocephala within Syndermata. Here, we present the first study of the phylogenetic affinities of Chaetognathifera with genomic representation from all major phyla, and explore the potential of macrosynteny to better understand these relationships. For this latter aspect, we also developed a new jackknifing procedure to assess the robustness of linkage groups inferred by macrosyntenic analyses. We show that the phylogenetic relationships between these clades are corroborated through a variety of gene selection and analysis methodologies. This provides clear evidence of Acanthocephala as a derived clade within Syndermata as sister to Seisonidea, and that Gnathostomulida is sister to Syndermata+Micrognathozoa, with Chaetognatha as the earliest diverging clade within Chaetognathifera. On the other hand, we found that macrosyntenic patterns cannot resolve this question. Moreover, almost all possible linkage groups involving chaetognathiferan species lack robusticity and hence, should not be considered reliable. As a consequence so far, in Chaetognathifera none of the bilaterian ancestral linkage groups can be reliably found and independent massive chromosomal rearrangements occurred. We therefore strongly suggest that studies of macrosynteny should not only assess the significance of possible linkage groups, but also the robusticity of these linkage group inferences. Furthermore, we also present a script for this purpose, which can be found at: https://github.com/JFFleming/MacrosyntenicJackknife

Cities are characterized by elevated temperatures, increased pollution, and high-density human populations which often are accompanied by changes in available resources, like food. These shifts have the potential to drive phenotypic divergence in urban wildlife. Functional morphological traits, like body size, can mediate interactions between wildlife and habitat and are closely tied to life history and fitness. While examples of functional morphological variation associated with urbanization are increasing, variation in such traits as a response to urbanization remains unexplored for most taxa. Here, we investigated morphological divergence between urban and rural populations of house mice (Mus musculus domesticus). House mice are globally distributed in diverse habitats and are a model system with a wealth of phenotypic data, making them useful for the study of the impacts of urbanization on morphology. Using a paired replicate design, we sampled urban and rural populations in three distinct metropolitan regions in the eastern United States. We found that body size was smaller in urban populations. Using 3D geometric morphometrics, we also analyzed variation in cranial shape across habitats. Differences in cranial shape were largely allometric, that is, driven by differences in body size. However, we also uncovered evidence of cranial shape variation between habitats not explained by size. In contrast, we did not find evidence for habitat-driven differences in cranial capacity independent of size. Overall, our results suggest a key role for body size in mediating morphological responses to urbanization and highlight the potential of house mice as a globally-distributed model for urbanization.

For over four decades, the bivalve Anomalocardia flexuosa has been recorded in Hong Kong coastal waters. However, the known native distribution of this heavily exploited commercial species is restricted to the Atlantic coast of South America, raising questions about the biogeographical validity of the Hong Kong populations. By employing an integrative taxonomic approach combining morphological re-evaluations and molecular phylogenetic analysis of the COI gene, we confirm that the species in Shui Hau, Hong Kong, China, has been historically misidentified. The population belongs to Cryptonema producta (syn. Anomalocardia producta).

In odonates, hybridization has been documented in several species pairs, yet the long-term persistence of hybrid phenotypes remains poorly understood, particularly in Neotropical systems. Here, we investigate a putative hybrid zone between Ischnura capreolus and Ischnura cyane, which overlap along an altitudinal gradient in the Colombian Andes. We combined temporal sampling across sympatric and allopatric localities with multilocus genetic analyses and geometric morphometrics of wings and male caudal appendages to evaluate patterns of admixture and phenotypic intermediacy. Morphometric analyses revealed pronounced differentiation between I. capreolus and I. cyane, especially in allopatric populations, whereas individuals from sympatric localities displayed increased morphological overlap and intermediate phenotypes. Genetic analyses based on nuclear and mitochondrial genes identified two main genetic clusters corresponding to the parental species, with evidence of admixture and shared haplotypes in sympatric localities. Patterns of differentiation varied among loci, with mitochondrial markers showing broader haplotype sharing than nuclear loci. Across sympatric localities, morphologically intermediate individuals remained consistently present through time, particularly in Anolaima, where they represented the dominant phenotype across sampling periods. These results support concordant morphological and genetic evidence consistent with persistent hybridization and introgression between I. capreolus and I. cyane in their zone of contact and suggest that incomplete reproductive isolation contributes to the long-term maintenance of admixed phenotypes in this hybrid zone.

The delineation of closely related species remains a persistent challenge in Zoantharia, where morphological plasticity and limited genetic differentiation complicate taxonomy. In this study, we investigated the phylogenetic relationship between the widely distributed sibling taxa Palythoa tuberculosa (Indo-Pacific) and Palythoa caribaeorum (Atlantic) using ultraconserved elements (UCEs) recovered from genome skimming. A dataset comprising 116 loci (35,699 bp) across 37 specimens from Brazil, the Red Sea, Okinawa, and New Caledonia was analysed using both concatenated maximum-likelihood and coalescent-based approaches. Phylogenetic reconstructions did not recover monophyletic relationships corresponding to either species or geographic origin, instead revealing intermixed lineages across the Indo-Pacific and Atlantic regions. Concordance factor analyses indicated low gene concordance and moderate site concordance, suggesting pervasive gene tree discordance rather than a lack of phylogenetic signal. These patterns are consistent with previous studies based on mitochondrial, nuclear, and reduced-representation datasets, indicating that increased marker resolution does not resolve species boundaries within this complex. The observed lack of differentiation may reflect ongoing or recent connectivity among populations, potentially facilitated by long-distance dispersal promoted by anthropogenic rafting or historical range expansion, biological invasion, or biological processes such as incomplete lineage sorting. The results support the hypothesis that P. tuberculosa and P. caribaeorum represent a species complex or a case of incipient speciation rather than fully distinct evolutionary lineages. These findings indicate that genome-scale data alone may be insufficient to resolve very recent divergences, supporting the need for integrative approaches to resolve complicated species boundaries in zoantharians.

True crabs (Brachyura) are among the most iconic marine arthropods, representing noteworthy examples of morphological and ecological disparity. A striking feature of brachyurans are their anterior pincer-like appendages: chelipeds. These structures showcase a large diversity of morphologies that reflect ecology and overall multifunctionality. Yet, a comprehensive assessment of appendage functional morphology within phylogenetic and ecological trait contexts has never been attempted. By combining 3D geometric morphometrics, finite element analyses, multilocus molecular phylogeny, and ecological trait data for 80 crab species, including three fossil forms, we unveil a complex evolutionary history for crab chelipeds. Despite extreme shape diversity amongst chelipeds, stress distributions are very similar across taxa and hint a many-to-one pattern. High concentrations of chelipeds within constrained morphospace regions associated with peak pinch forces illustrates that brachyuran morphologies optimised for shell crushing may have arisen in the Cretaceous. Deviations from this morphospace highlight the diversification of non-shell-crushing life modes and the influence of sexual selection on appendages. Neither cheliped shape nor pinch force show phylogenetic signal. Together these results indicate that the evolution of cheliped shape is closely associated with, and inferred to have been strongly influenced by, crab ecology, biomechanical needs and sexual selection. SIGNIFICANCE STATEMENTChelipeds, the pincer-like claws of crabs, are among the most morphologically diverse appendages within Arthropoda, yet the evolutionary forces driving this diversity remain poorly understood. By integrating 3D geometric morphometrics, biomechanical modelling, molecular phylogeny, and ecological data across 80 crab species including fossil forms, we demonstrate that cheliped morphology is driven by ecology, biomechanical demands, and sexual selection rather than phylogenetic relatedness. The multifunctionality of these structures produces strong evidence for many-to-one mapping of form to function. Morphologies optimised for durophagy appear to have originated in the Cretaceous, with subsequent diversification into manipulative and sexually selected forms from a morphologically flexible foundation. These findings demonstrate that cheliped diversity reflects a complex interplay between ecological specialisation, biomechanical optimisation, and sexual selection across Brachyura.